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Abstract 

The effectiveness of viscous elements in introducing damping in a structure is a function of 
several variables, including their number, their location in the structure, and their physical 
properties. In this paper, the optimal damper placement and tuning problem is posed to op- 
timize these variables. Both discrete and continuous optimization problems are formulated 
and solved, corresponding, respectively, to the problems of placement of passive elements 
and to the tuning of their parameters. The paper particularly emphasizes the critical com- 
putational issues resulting from the optimization formulations. Numerical results involving 
a lightly damped testbed structure are presented. 

1. Introduction 

A problem of considerable importance in the development of technology for future space 
structures is the analysis and optimization of passive elements placed in these structures. 
Passive damping introduced by these devices is an effective mechanism for reducing peak 
responses in the vicinity of resonant frequencies for lightly damped systems. This not only 
enhances the stability of the open-loop system, but also allows for the implementation of 
more aggressive control strategies to achieve greater performance. This philosophy is being 
pursued on a series of Control Structure Interaction (CSI) testbeds at the Jet Propulsion 
Laboratory. 

The effectiveness of viscous elements in introducing damping is a function of several vari- 
ables, including their number, their location in the structure, and their physical parame- 
ters, namely damping and stiffness coefficients. This paper is concerned with the optimal 
placement and tuning problem for the passive viscous dampers with emphasis on its com- 
putational aspects. 

Two qualitatively different optimization problems arise in this context: a combinatorial 
optimization problem which determines the placement of elements, and a mathematical 
programming problem which optimizes (tunes) the damper parameters. In our approach, a 
simulated annealing strategy [4] is used for the combinatorial optimization problem, while a 
sequential quadratic programming algorithm (SQP) [2] is applied to the damper parameter 
optimization problem. One of the most important ingredients in any optimization prob- 
lem is the cost functional evaluation, regardless of the performance metric that is used. 
This is particularly true for the optimal damper placement and tuning problem due to the 
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complexity of the system. The performance metric chosen here is the W 2 -norm of selected 
transfer functions of interest. An excellent candidate is the transfer matrix between external 
disturbance inputs and the controlled outputs. 

It is well known that the computation of the 'W 2 -norm requires solving a Lyapunov equation. 
However, due to the high-dimensionality of the system model, it is unrealistic to use the 
full-order model in any computation. A reduced-order model must be generated to make the 
computation involved more manageable. The Ritz reduction method that has been studied 
in [1] is employed to reduce the numerical bottleneck created by solving large systems of 
this type. 

The paper is organized as follows. Section 2 presents the dynamic model of a viscously 
damped structure. The general optimal damper placement and tuning problem is formu- 
lated in Section 3 with a review on the computation of the W 2 *norm of the particular transfer 
matrix which is chosen to be our performance metric. Section 4 addresses the computation 
issues involved in our optimization problem. In particular, the Ritz reduction method will 
be described in detail. A number of numerical examples involving the JPL testbed structure 
are presented in Section 5. Finally, concluding remarks on future work are given in Section 
6 . 


2, Dynamic Modeling for Viscously Damped Structures 

Throughout this paper, it is assumed that the dynamics of the undamped structures can 
be described by a linear, second-order matrix differential equation of the form: 


Mz + I\ z = Bdd . 


(1) 


Here z denotes the n-dimensional vector of generalized coordinates, d is an /-dimensional 
external forcing input vector, M is the n x n symmetric, positive definite mass matrix, K 
is the n x n symmetric, positive definite stiffness matrix, and Bd is the n x / forcing input 
influence matrix. 

Assume that a discrete passive damper is placed between two nodal points in the structure, 
replacing the original structural element. The passive damper is modelled as a device that 
applies a force at the nodal points with equal magnitude but in opposite directions and 
proportional to the relative displacement and velocity between the nodal points. 

The dynamic structural model incorporating the damper actuator force, u, is written as 

M2+ Kz = bu + B d d (2) 

where the vector b represents the influence vector associated with u. The force u generated 
by the damper is modelled as a constant linear combination of collocated position and 
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velocity feedback so that 


u = ~(k p y p + k v y v ) (3) 

with y p = b T z and y v = b T z where y p and y v denote the position and velocity “mea- 
surements,” respectively, and k v denotes the damping rate, which is always taken as a 
nonnegative quantity to ensure stability. The parameter k p is only required to be greater 
than or equal to the value -k e , where k e denotes the stiffness of the structural element that 
has been replaced by the damper. When -k e < k p < 0, the structure is softened, while 
k p > 0 causes the structure to be stiffened. 

Hence, the dynamic structural model with the inclusion of a passive damper can be repre- 
sented as 

Mz + (K + k p bb T )z = b( -k p b T z - k v b T z) + B d d , (4) 

or 

Mz + (k v bb T )z + (K + k p bb T )z = B d <l . (5) 

A more general model including multiple passive dampers can be written as 

Mz + (£ k v b,bi T )* + (A' + £ k p ,bibi T )* = B * d (6) 

*=1 1=1 

where n p is the number of passive dampers in the structure. 

3. Optimal Placement and Tuning Problem for Passive Dampers 

The general optimal placement /tuning problem of passive dampers can be posed as 

min min J cost {B p , K py K v ) (7) 

K p eK Vt K v e£v B p £Bp 


where 


• Jco»t{Bp,K p ,K v ) is defined as the performance metric for the optimization with a 
given damper configuration of locations corresponding to B p and the corresponding 
stiffness and damping rate I\ P and K v 


• Bj> — {(frj, i ^137 • *1 » *2* • • • 1 *n r € Afv > ia 7^ ’ Vo , / 3 1,2,. 

( b ia is the influence vector corresponding to the i Q th location). 


--,n p (a / /?)} 


JC-p={(k Pii ,k Pij ,...,k Pirt ') : e £ */ f,Va,/?= 1,2 ,...,n p (a ± /?)} 

( k Pj is the stiffness correction corresponding to the damper at j th location, and 


k. ^ kn -f- k e ^ k 3 

3 mtn — P} ' e J — 

where k e> is the element stiffness of the undamped structure at j th location, k smin and 
ksmax are the lower and upper bound of the damper stiffness). 
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• )C v ={(k V ' i ,k v ,' i ,...,k Vin ') : i n , € Afp, i a / ^,Va,/3 = l,2,...,n p (a / /?)} 

(/t Vj is the damping rate corresponding to the damper at j th location, and 

0 5: k V j — k VmaI 

where k VmmM is the highest possible damping rate for the passive damper). 

* J\f-p is defined as the set of all candidate locations for placement. 


It is clear that the above optimization problem is a joint “continuous+discrete” optimiza- 
tion problem. The selection of locations ( B p ) for placement is a “discrete” combinatorial 
optimization problem while the selection of values for A p and A„ (tuning) is a continuous 
mathematical programming problem. 

Two types of performance metrics are typically considered. The first one is the structural 
modal damping for selected modes. The computation involved is to solve for the eigenvalues 
of the “A”-matrix obtained from writing (6) in first-order form for a given damper config- 
uration with corresponding damper stiffness and damping coefficients. The second type 
of criterion requires both the external disturbance input vector and the controlled output 
vector to be specified. As discussed in [6], a meaningful and numerically tractable criterion 
for the associated optimization problem is to minimize the W 2 -norm of the transfer function 
from d to y a . In addition, a weighting function W'd(s) can be used to model the spectral 
property of d and a weighting function W p (s) can be used to improve the performance of 
y„ over a certain frequency range. In this case, the cost functional is simply 

Jco,i = ||Wp(s)G>(s; B„,K P , A’„)W',(5)|| 2 ( 8 ) 

where G p (s; B p , A' p , K v ) is defined as the transfer matrix from the d to y 0 with a given 
damper configuration of locations corresponding to B p and with the corresponding stiffness 
and damping coefficients, A' p and K v . For a given damper configuration (B p , A' p , A'„) and 
the weighting functions (IT p (s), IT d (s)), the 7f r norm can be computed through the solution 
of a specific Lyapunov equation. 

Define 

T(s) = W p (s)G p (s; B p , K p , K v )W d (s) 

and assume that T(s) has the state-space realization (A, A, C ) where the matrix A is asymp- 
totically stable. Then the corresponding 7f 2 -norm of T(s) is simply 

||T’(s)|| 2 = [trace(CPC T )} 1/2 = [trace(B T QB )] 1/2 

where P and Q are the positive semi-definite solutions of the following two Lyapunov 
equations: 

AP + PA t + BB T = 0 (9) 
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and 


A t Q + QA + C T C = 0 , 


( 10 ) 


respectively [3]. 


4. Computational Issues and Model Reduction 

As discussed in the previous section, the damper placement and tuning problem includes 
solving a nonlinear mathematical programming problem for tuning, and a combinatorial 
optimization problem for placement. 

In particular, the combinatorial optimization problem is known to be difficult due to the 
fact that the potential number of candidate locations for placement (N p ) will be large in 
large space flexible structures. However, relatively few passive devices (n p ) will be available. 
In general, N p >> n p , and the total number of combinations, , (Ar'l w ,)? ’ usually very 
large. Therefore, it is impractical, if not completely impossible, to try the exhaustive search. 

In our approach, a sequential quadratic programming algorithm (SQP) [2] is applied to 
the damper parameter tuning problem while a simulated annealing strategy [4] is used for 
the combinatorial optimization problem. The question of developing a hybrid approach for 
combining these strategies into a single approach will not be dealt with here and is one of 
our future research topics. 

Our current strategy is to solve each of these problems individually. One approach is to solve 
the damper parameter tuning problem for each candidate location first. These parameters 
will then be used to evaluate the cost functional in the simulated annealing process. 

Another approach is to use a “pruning” process after each of the candidate locations is 
“tuned.” This pruning process is simply to choose the top N p candidate locations accord- 
ing to the ranking of their respective optimized cost functional where N p >> N p > n p . 
An exhaustive combinatorial search is then conducted throughout this subset to find the 
“optimal” combination of elements which yields the smallest 7^ 2 -norm cost. This ad hoc 
pruning approach has been demonstrated to be quite useful. However, it is difficult to make 
a general statement regarding the solutions of these sub-optimal approaches as compared 
to the optimal ones. 

As stated in the Introduction, one of the most important ingredients in any optimization 
problem is the cost functional evaluation. This is particularly true for the optimal damper 
placement and tuning problem due to the complexity of the system. The performance 
metric chosen here is the W 2 - norm of selected transfer functions of interest. 

The procedure to compute the H 2 *norm of a stable transfer matrix has been given in Section 
3 and requires solving a Lyapunov equation. However, it is impractical, if not impossible, 
to use the full-order model in the computation of the W 3 -norm since the order of the model, 
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2 x n, is typically very large. Hence, a high-fidelity, low-order, reduced model must be used 
to perform the required computation efficiently. 

The Ritz reduction method that has been studied in [1] is employed to reduce the numerical 
bottleneck created by solving large systems of this type. Details of this model-reduction 
method will be described in the rest of this section. 


The Ritz Reduction Method 


To solve the optimization problem posed in the previous section, it is impractical, if not 
impossible, to use the full-order model in the optimization process since the order of the 
model, 2 x n, is typically very large. Hence, a high-fidelity, low-order, reduced model must 
be used to perform the required computation efficiently. 


The model-reduction method considered here is a second order reduction technique based 
on reducing the number of generalized coordinates of the system via a transformation of 
the form z — Vq, where q £ R N with N < n. Applying the transformation V to (6) results 
in the reduced-order model 


C V T MV)q + 


'£k v ,(V T b 1 )(V T b t ) T 


«=i 


<z + 


(v t kv) + X>,.( r T b i )(v T b i ) 1 


* = l 


q = ( V T B i )d . 

( 11 ) 


The transformation matrix, V, consists of the first m (m << n) eigenvectors corresponding 

to the first m eigenvalues, {uq, u? u; m }, and an additional Ritz vector to account for 

the static correction for each of the forcing inputs. This method will be referred to as the 
“Ritz reduction method.” A detailed discussion on this subject can be found in [1]. 

Suppose that the lowest m eigenvalues and their corresponding eigenvectors are known 
and 4> m is defined as the n X m matrix consisting of the m eigenvectors corresponding to 

{uq, w 2 , t^m}- Then the desired Ritz vector corresponding to b t ( i th damper) is simply 

the solution to the following linear equation: 


K^i = b, . 


It is desirable for the transformation matrix to preserve A/ orthonormality. Therefore, 
needs to be M-orthonormalized. This is done easily by first 


1. making ipi M-orthogonalized to 

= x pi - $ m ($ m T Mxl>i) 


( 12 ) 


and then 
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2. making M-normalized, i.e., 


(13) 


<t>r = • 

Similarly, the desired Ritz vector corresponding to b dj (j th external disturbance input) can 
be computed using the same procedure. 

Note that for each of the forcing inputs, one Ritz vector needs to be computed. The 
forcing inputs could be either the force inputs corresponding to the dampers or external 
disturbance inputs. Let <£• denote the M-orthonormalized Ritz vector corresponding to 
the influencing input vector, 6,-, and <i> r di denote the M-orthonormalized Ritz vector 
corresponding to the j th external disturbance influencing input vector, b dj . Note that each 
of the corresponding Ritz vectors is M-orthonormalized to 4> m ; however, the (n p + /) Ritz 
vectors may not be M-orthogonal among themselves. An additional M-orthogonalized step 
is required. Define 

= [4>\ & - 4> r n, 4>\4> T ^ -<&,] 

and form 

M r itz = and K riu - $* iu K$ riu 

to find 4> r<(2 such that $ rilz is M ri( ,-orthonormalized, i.e., 

= I(n p +l)x(n r +l) i and A r i«z ^riti 

where fl r ,n = diag[u> ri u; r3 ... w rnf+1 ]. 

Define 4> ritz = $ ri <z * 4> r ii*, then the M-orthonormal transformation matrix V is 

V = \$ m $r,u] 


and Eq. (11) is equivalent to 

InxnQ + XB T bi){V T bi) T )q + 


i=l 


ni + ^k p .(v T b,)(v T b i y 


i=i 


q = ( V T B d )d (14) 


where N — m + n p + l is the order of the reduced model, and fijv = diag[Q. m Q ritz ]. 

The reduced-order model in Eq. (10) can also be rewritten in the state-space representation 


as 


x = 


0 JVxN -LvxN 

-& N - k P xv T b,){V T b i ) T -YZi a KX^b^bd 1 


X + 


0 JVxn 

V T B d \ 


d (15) 


where x = 


is the state vector. 
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5. Numerical Examples 


A detailed description of the JPL testbed can be found in [5] (see Figure 1). Briefly, the 
system is modeled with 249 degrees of freedom and contains 186 candidate locations to 
insert passive damping elements. 

Because the accuracy of the cost functional evaluation methods is of paramount importance 
in the optimization process, Table 1 contains a comparison of eigenvalue approximations 
using the full-order model, the Ritz reduced model, and a modally reduced model. The 
second column in part (a) of the table contains the eigenvalues of the undamped nominal 
system. All of the other values correspond to the damped system with three viscous dampers 
placed at the locations 132, 140, and 142. It is assumed that the three dampers have the 
same damping and stiffness coefficients: 320 lbs - sec l in and 8,000 lbs /in respectively. 

The conclusion here is that the Ritz reduction method yields high-precision estimates with 
enormous reduction in computation. In this example, instead of solving a 498 x 498 eigen- 
value problem, the results can be obtained by solving a 30 X 30 eigenvalue problem which 
results from the Ritz reduction method. However, the modally reduced model produces 
inaccurate results. What is of equal significance is that not only does the modally reduced 
model produce inaccurate results, it also leads to inaccurate trends for choosing damper pa- 
rameters. Figure 2 contains damping predictions of the second system mode as a function of 
the damper viscous parameter coefficient. Note that the full and Ritz reduced models lead 
to an optimal coefficient of approximately 500 lbs — sec/ in, while the modally reduced model 
leads to a significantly larger value that is far from optimal. The Ritz reduction method 
also leads to very accurate approximation to the 7i2-norm, with 6 digits of accuracy. 

Table 2 contains the eigenvalues of the damped system where the three dampers are placed 
at the locations 6, 19 and 91. The three locations are the simulated annealing solution to 
the optimal damper placement problem. The performance metric is the W 2 -norm of the 
transfer function from an input disturbance located at grid point 412 between the third 
and fourth bays of the structure, to the outputs consisting of all of the nodal displacements 
directly beneath the trolley (see Fig. 1). The disturbance was generated as the output of 
a 6 <h -order low-pass filter with a bandwidth of 25 Hz. This weighting function is chosen 
to reflect the objective of disturbance reduction in the frequency range below 25 Hz. A 
representative comparison of the undamped and damped frequency responses is given in 
Figure 3. 


6. Concluding Remarks 

The use of strategically placed and tuned passive elements in future large space structures 
will play a significant role in their design and development. The ability to analyze, predict, 
and ultimately optimize system performance with respect to these passive devices is critical 
for the application of this damper placement technology. 



A comprehensive overview of the optimal damper placement and tuning problem was pre- 
sented in this paper. Approaches and computational aspects of the associated optimization 
problems were discussed. The results of the paper indicate that significant levels of damping 
can be introduced into these structures in a very systematic and tailored fashion. 

Although reasonably good results have been demonstrated using the approach presented 
here, the combined discrete plus continuous optimization problem was essentially solved for 
each individually. This is the major drawback of our current approach. Our future work will 
concentrate on the development of a hybrid approach to jointly solve the two qualitatively 
distinct optimization problems. 
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Figure 1 . JPL CSI Phase B Testbed 
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Mode 

Undamped System 

Damped System 

249 Modes 
( Full order) 

12 Modes plus 
3 Ritz vectors 

15 Modes 
(Truncation) 

1 





2 

5.4263 

5.2940 


5.3262 

3 

7.4565 



6.9540 

4 

11.6777 

10.4862 

10.4862 

10.4493 

5 

17.4248 

17.4386 

17.4386 

17.3444 

6 

20.8423 

20.8236 

20.8236 

20.7055 

7 

31.1387 

31.2231 

31.2231 

31.0481 


(a) Frequency (in Hertz) 


Mode 

Damped System 

249 Modes 
(Full order) 

12 Modes plus 
3 Ritz vectors 

15 Modes 
(Truncation) 

i 

0.0179 



2 

4.5744 

4.5744 


3 

25.5358 

25.5357 


4 

32.6380 

32.6379 


5 

0.9033 

0.9034 


6 

1.3197 

1.3197 


7 

0.5013 

0.5016 

0.5031 


(b) Damping (in %) 


Table 1. Undamped and Damped Eigenvalues 
(Damper Locations: 132, 140, and 142) 


Mode 

Frequency (Hz) 

Damping (%) 

i 

IHHB5QEH 

0.0245 

2 

■ 

6.8905 

3 

7.1748 

10.0192 

4 

11.4717 1 

3.7751 

5 

17.5924 

3.2823 

6 

20.9413 


7 

31.1573 



Table 2. Eigenvalues of the Damped System with Tio Optimized 
Damper Locations at 6, 19, and 91. 
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Figure 2. Damping Prediction by Reduction Methods 



Figure 3. Disturbance Frequency Responses of Undamped and Damped Systems 
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